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procedure is used which makes the kinetic equations employed exact. 
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1. Introduction 

In this letter, we should like to report on a result which can be obtained, for 
the model in question, without approximations from standard quantum theory of 
open systems governed by the linear Liouville - von Neumann equation. Irrespec- 
tive of this, it contradicts the second law of thermodynamics. Leaving details of 
physical motivation to a next publication, we mention here just the fact that the 
main inspiration for construction of the model has been taken from biology, namely 
from topological changes of biologically important molecules upon detecting, at a 
specific site (receptor), particles (excitations, molecules or molecular groups) to be 
processed [|TJ. 

Previous version of the model has been published in ||. Detailed treatment of 
this model || as well as other microscopic models of open quantum systems work- 
ing on analogous principles ||, || [F[ |] revealed property of spontaneous (i.e. not 
induced by external flows) self organization. This then leads to such unexpected 
phenomena contradicting basic principles of statistical thermodynamics as, e.g., vi- 
olation of consequences of the detailed balance (in connection with impossibility of 
rigorous justification thereof). From this, implicit violations of the second law of 
thermodynamics could be deduced as mentioned in, e.g., || ||. Here, we recon- 
struct the model so that it is able to work cyclically and without compensation as 
a perpetuum mobile of the second kind, in the sense contradicting explicitly the 
Thomson formulation of the second law of thermodynamics ||. 

2. Model 

The fully quantum Hamiltonian of our model in its simplest version (see also 
[T0| ) can be as usual written as a sum of the Hamiltonians of the (extended) system, 



thermodynamic bath, and that of the system-bath interaction. Thus, 

H = H S + H B + H S _ B (1) 

where 

H s = ^(c^co + CqC_i) ■ \d)(d\ + I{c\c Q + cjci) ■ \u)(u\ + K(c\c^i + cL±ci) 

+ € -[l-2clc }-[\u}( U \-\d}(d\}. (2) 

As for the Hamiltonians of the bath Hb and that of the system-bath interaction 
Hs-b-i they are not important here. Before explaining the symbols used in (Q), 
let us only add that we shall assume H$~b to consist of two additive and non- 
interfering contributions H' S _ B and Hg_ B . Here H' S _ B causes transitions between 
states of the central system \u) and \d) (see below) while Hg_ B is responsible for 
a sufficiently intense (and, for the sake of simplicity, equally strong) dephasing of 
the particle states at all three sites. As an example, one can take a model of 
the bath consisting of harmonic oscillators Hb = J2k hookb\bk interacting with the 
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central system by a linear site- local coupling causing the \u) <->• \d) transitions 
exactly as in || but complemented, for our purposes, by, e.g., a non-interfering 
site-local coupling of the particle to the bath with the same coupling constant at 
all three sites. This means, e.g., H S _ B = -^J2k^kG k (b k + bl_ k )[\u)(d\ + \d)(u\] 

+ jf Ekk' hy/WkWk'9kk'(h + b ] _ k ){b k , + ftLfc/) (cliC-i + cjc + c\ x c +1 ) = H' S _ B + H%_ B . 
These specific forms of H B and Hs~b will not be, however, used below. 

As for H s in ([|), we have chosen the simplest model with only three particle 
states, i.e. sites (m = —1,0 and +1). Operators and c m , m = —1,0, or +1 des- 
ignate the particle creation and annihilation operators at site m. For simplicity, we 
shall assume only one particle in the system. This is why we do not need the com- 
mutational (anticommutational) relations of the particle creation and annihilation 
operators. In (H), I, J and K are transfer (hopping or resonance) integrals connect- 
ing the sites involved. Worth noticing is that in ([|), the forth and back transfers 
in any pair of the sites are always with the same amplitude as a consequence of the 
hermicity of Hs- The one-directional character of the process reported here is not 
owing to a contingent difference between these amplitudes but results exclusively (as 
it will become clear later on) from the existence of spontaneous processes. Site '0' is 
understood to be attached to a central system representing, e.g., a specific molecule 
or molecular group (e.g. tail connecting site '0' with either site '-1' or '+1' but not 
both simultaneously). This central system is assumed to have (in a given range of 
energies of interest) two levels with energies ±e/2 (with corresponding states \d) 
and | it); we assume e > 0). At the moment when the particle is transferred to site 
'0' attached to the central system, the relative order (on the energy axis) of the two 
levels of the central system gets interchanged. (This of course causes instability of 
the central system with respect to the \d) — > \u) transition which is the same effect 
as additional-load-induced instability of a ship in water.) Asymmetry in transfer 
rates \u) <-» \d) is ensured by the above spontaneous processes with respect to the 
bath. Conversion of this up-and-down asymmetry into the left-and-right one is then 
owing to the above special form of the first two terms on the right hand side of (0) 
discussed below. Technically, this imbalancing is due to the fourth term on the right 
hand side of fl2|) proportional to 1 — 2cqC and may be in reality due to correlation 
effects as the particle transferred may, upon its transition to site 0, change the topol- 
ogy (originally stable conformation may become energetically disadvantageous) or 
orientation of the central system in space. (Such changes condition activity of many 
biologically important molecules in living organisms |I|].) 

3. Equations of motion 

In our theory, we shall closely follow || but shall be interested in just a station- 
ary situation. That is why we can apply any proper kinetic method yielding time 
derivatives of matrix elements of the density matrix of our extended system (particle 
+ the central two-level system), and to set these time derivatives zero. In order to 
avoid unnecessary technical complications, we have avoided time-convolution meth- 
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ods. From mutually equivalent time-convolutionless methods, we have chosen (like 
in ||) that of Tokuyama and Mori |Tl|j . In order to avoid unnecessary discussions 
about role of approximations, we have employed one of the scaling procedures turn- 
ing, in the limiting sense, Tokuyama-Mori equations into exact kinetic equations. 
This means the following steps: 

• We formally introduce a joint small parameter, say g, of both H$-b and all 
the hopping integrals J, J, and K in the sense of setting Hs-b oc g but 
/, J, K oc g 2 . Right here, let us remind the reader of the fact that this scaling 
of some parameters from Hs is what distinguishes our approach from scaling 
standardly used in weak coupling theories (inapplicable to our presumably 
intermediate or rather strong coupling situation discussed here) |12| . 



• We introduce new time unit r oc g 2 , i.e. introduce new time t' = t/r . This 
step formally disappears as far as we are, like here, interested in just the 
stationary (long-time) asymptotics. 

• We divide our general Tokuyama-Mori kinetic equations by g 2 and perform 
the limit g — > 0. 

The reader can easily see that this method preserves just the lowest order (in g 2 ) 
terms which can be calculated exactly. The result can be reported as follows (for 
details see a next extended publication): 

First, we arrange all the 36 matrix elements p ai (t) of the density matrix of our 
extended system (a, (3... = md or mu with m = or ±1 while d or u designate the 
states of the central system) in groups of nine designating 



P-\u-\u-> POuflu, Plu,lu, P-luflu) P-lu,lu) POu-lui p0u,lui plu,~lu, Pluflu I , 



{Pddf 



ud I 



P-ld,~ld, POdfld, Pld,ldi P-ld,0d, P-ld,ldi POd-ld, P0d,ldi Pld,-ld, Pldfld J , 
P-lu-ld, POufid, Plu,ldi P-lufldi P-lu,ld, P0u-ld, P0u,ld) P±u-ld, Plufld 



(3) 

and similarly for pd u - Superscript • • - T designates transposition. Then the above 
kinetic equations in the asymptotic time domain (and after the above scaling) read 
as 
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(4) 
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Here, in the square matrix, all the elements are in fact blocks 9x9. Hence, the whole 
set splits into two independent sets of 18 equations each; we shall be interested just 
in that one for p uu and p^d- This reads as in (^) with typical forms of the block 9x9 
matrices A, B, C and V. Here 

A = 



( 
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(5) 

Here k = —ie/h — 0.5(1^ + Tj) and • ■ •* designates complex conjugation. and 
designate bath-assisted uphill and downhill transfer rates in our two-level system 
calculated as if, formally (owing to the above scaling) I = J = K = 0. Let us 
mention that via I\ and Tj, the (initial) bath temperature T enters the game. 
These rates are known from the usual Pauli Master Equation approach to general 
kinetic problems and, consequently, are related to each other by the detailed balance 
condition 1^ = I^e - ^, (3 = 1/(&bT). In general, however, owing to intermixture of 
coherent and incoherent transfer channels, the detailed balance condition does not 
apply to our particle transfer problem. Finally, 2T is the dephasing rate determined 
by Hg_ B . As for the block B, it is fully diagonal with diagonal elements Bn, . . . , £> 99 
equal to I>, T u T h (I> + TJ/2, T h (I> + 1^/2, (I> + TJ/2, T T and (I> + TJ/2, 
respectively. Next, 

V = 
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(6) 

As for the block C, it reads as B except for the interchange I\ <-> The form of 
all the matrices in (fD is exactly the same as, e.g., that one which we would get from 
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the stochastic Liouville equation SLE [|13j provided, however, that Hs-b is replaced 
by a proper stochastic (e.g. Gaussian delta-correlated) potential field acting on the 
central system. The only difference between our form of the A — T> blocks and 
that of the same matrices in the corresponding SLE theory is dictated by physics of 
the problem: Namely, in contrast to the SLE approach, spontaneous processes with 
respect to the bath naturally appear in our fully quantum model (Hs-b is, in our 
case, a coupling to a genuine quantum bath). Thus, < or even <C in our case. 
As SLE is quite standard and well understood, this comment will hopefully turn 
attention of suspicious readers from potential speculations about correctness of our 
approach to the form of our Hamiltonian. It is the latter what is responsible for the 
striking conclusions obtained. 

4. Stationary flow and power output. 

Upper half of (Q) provides a set of 18 linear algebraic homogeneous equations 
for 18 components of the density matrix p uu and p^d (remaining components be- 
coming decoupled and irrelevant in what follows) of the rank 17. So, it should be 
complemented by the normalization condition 

+i 

^ t [Pmu,mu Pmd.md\ L (7) 
m=— 1 

Then this set and (^) determine the relevant components of the density matrix 
uniquely. From the physical meaning of the transfer rates, one can then determine 
the flow '-lWOW+l'^'-l' as 

J — ^lP0d,0d — TiPom.Om- (8) 

Calculated flow is, as simple inspection or numerical solution shows, always nonzero. 
Already this implies, owing to the persistent character of the flow, important conse- 
quences. We are, however, here interested rather in the possibility of converting the 
heat from the bath into a usable work. For that, we connect the particle running 
through our system organized as a circle of three sites (as '-1— >0— H-l— with a 
screw (with its axis perpendicular to the plane of our three sites '-1', '0' and l +V). 
The above persistent flow could then drive the screw at the cost of just the thermal 
energy of the bath, converting thus the latter directly to the mechanical work. So 
the existence or nonexistence of the stationary flow according to the above defini- 
tions converts to a question of a direct violation of the 2nd law of thermodynamics 
as applied to our system in the sense of the existence or nonexistence of the per- 
petuum mobile of the second kind. The problem can be easily solved numerically in 
our model. The idea is to provide the system with potential steps the particle (and 
the screw) must overcome when passing from '-1' to '0', from '0' to '+1' or from '+1' 
to '-1'. (When passing virtually in the opposite directions, the potential steps as felt 
by the particle would get the opposite sign.) We have chosen these steps equal, des- 
ignating their value as A/3. Hence A > is the mechanical work the particle exerts 
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Figure 1: Thermal-to- mechanical energy conversion power-output W in units I 2 /h 
(as determined by fllUp) as a function of A; / = J > 0, e = 50/, I\ = 0, Tj = 101 /h, 
2F = 3I/h. Parameter K = 0.1/, 0.05/ and 0.01/ for curves a, b, and c, respectively. 

on the screw during one turn in the above direction. One should realize that taking 
these steps as those of the potential energy of the screw connected with the particle, 
the latter potential energy is not unique as a function of the particle position on the 
above triangle of the sites. Instead, it is unique as a function of the rotation angle 
(p of the particle (or screw). In other words, when the particle performs one turn 
-1'— >'0'— 0^0 + 2tt, the potential ascribed to the particle connected 
with the screw increases by A though the particle formally returns to the same site. 

Technically, inclusion of the above potential steps is simple. As it follows from the 
above formalism or already even from the Liouville equation for the whole complex 
of the system and the bath, it only means to add to block matrices A and T> above 
the 9x9 block £ with all matrix elements equal to zero except for £ 44 = £ 77 = £ S8 = 
-£55 = -^66 = -£99 = ^A./(37i,). 

From the solution to 




A + £ B W Pun \ 

C V + £ )' [ Pdd J' 



(9) 

and (0), one can then determine the mechanical energy power output (as measured 
on the screw), i.e. the thermal-to-mechanical energy conversion power output (per 
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second) 

W = J-A. (10) 

Fig.l then unambiguously illustrates the positive result of our test whether the 
conversion of the thermal energy of our single bath to the mechanical work is, in 
our model, possible or not. One should realize the fact that the mechanical output 
is, as Fig. 1 shows, really positive. Hence, the system produces, owing to its 
selforganizational properties and properly timed opening and closing of the particle 
transfer channel across site '0', positive mechanical work. This work can, because 
of the construction of the model, be at the cost of just heat of the thermodynamic 
bath. Hence, we have a direct 'heat — > mechanical work' conversion. As we have just 
one bath, there is no compensation possible. Thus, our system working cyclically 
provides an example of perpetuum mobile of the second kind, violating thus the 
second law of thermodynamics in its Thomson form. 
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